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TECHNICAL NOTE NO . 1522 

CALCULATION OF UNCOUPLED MODES AND FREQUENCIES 
IN BENDING OR TORSION OF NONUNIFORM BEAMS 
By John C. Houbolt and Roger A. Anderson 

SUMMARY 


A procedure is presented for the calculation of frequencies and 
modes of nonunifoim Beams in uncoupled Bending and torsional vibration. 
Based on the principle of the Stodola method, the procedure oonsists 
of -solving the differential equation of equilibrium for vibration By 
a method of successive approximations. Basic principles of engineering 
Beam theory are employed in the method, and the integrations involved 
are perfoimed By improved numerical methods. An effort has Been made 
to perform all calculations in a manner consistent with the accuracy 
to which physiaal' constants in Built-up Beams are ordinarily known. 

Higher modes are readily found By use of the orthogonality 
relation Between normal modes. The frequency is found simply as the 
square root of the proportionality faotor existing Between modal 
deflection curves in successive approximations. All computations are 
tabular in form and are perfoimed mentally or with the aid of a slide 
rule. Comparison made with available oxaot analytical solutions shows 
that the method gives for practical purposes the. exaot answer. Special 
consideration has Been given to the treatment of various Boundary 
conditions that are found in the vibration of aircraft structures. 

The oantilever Beam, the free— free Beam, Beams with concentrated 
masses. Beams mounted on springs. Beams elastically ooupled to masses, 
and so forth, are shown to Be handled with practical simplicity. In 
order to servo as a guide in the solution of practical problems, the 
procedures for handling a number of different oases are illustrated By 
a liberal use of examples. 


INTRODUCTION 


In the dy n a m ic analysis of aircraft structures, the determination 
of the natural modes and frequencies is of Basio importance. A 
number of methods for calculating modes and frequencies have Been 
developed^ each method has certain desirable features. The objective 
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of the present paper is to develop a procedure which is readily 
learned hy anyone familiar with engineering "beam theozy, is easy and 
quick to apply without the use of complicated computing devices, and 
gives results within the range of accuracy with which the physical 
properties of the structure (mass, stiffness distribution, etc.) can 
ordinarily he determined/ An adaptation of -the successive- 
approximation procedure of Stodola (see reference l) fulfills these 
requirements and has the additional advantage that data necessary for 
the analysis of stresses due to vibration are obtained during the 
computations for the modes. In the present paper, this successive- 
approximation or iteration method 'is employed to obtain solutions to 
the differential equations of equilibrium for bending end torsional 
vibrations. 

The Stodola method, out lined, in reference 1, was originally a 
graphical-integration procedure for determining fundamental modes and 
frequencies. Burgees presented in reference 2 a numerical procedure 
for finding the fundamental frequency of a cantilever, which was in 
essence a series of approximate numerical integrations to determine 
the modal deflection followed by an energy solution for the frequency. 
Boukidis and Ruggiero (referrence 3) gave an application of the Stodola 
method, iri the form of a numerical procedure, which permitted cal- 
culation of the higher as well as the fundamental modes and frequencies 
of a free -free beam in symmetrical vibration. Inherent disadvantages 
in the 'method of higher-mode determination, as presented in reference 3 
however, 'have been pointed out by Beskin and Rosenberg in reference 4. 
These authors made Use of- the orthogonality relation between normal 
modes of vibration, in conjunction with Burgess's numerical procedure, 
to determine the higher modes of a cantilever. . . 

Two improved methods of numerical integration are presented in 
the present" paper for the solution of differential aquations by 
iteration. One of the methods is a summation process, somewhat ^ 
similar to the numerical procedure of reference 2; the other method 
■ is described by Newmark (reference 5 ) , who used it in beam and column 
■analysis and "called attention to its applicability 'to vibration 
problems. Both methods are simple to apply and. lead to relatively 
accurate modal deflections. The frequency is finally found to be the 
square root of the proportionality factor existing between modal- 
deflection curves in successive iterations. 

In order to perform the integrations, the boundary conditions 
on the vibrating member must be taken into account. Consideration 
has been given- herein to a number of. different boundary conditions. 

Both the symmetrical and the anti symmetrical modes of a free-free beam 
are treated in detail. Boundary conditions for the vibration of beams 
supported on springs, 'beams elastically coupled to masses, arid so 
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forth will he shown to he handled with -praotioal simplicity. Although 
the method can he applied to problems in which betiding and torsion are 
ooupled, the analysis will not he presented herein. 

70 , 

In order to determine modes higlier than the fundamental hy 
iteration, components of all modes lower than the. one to he determined 
must he removed. ' This operation is performed hy use of the orthogo— ' 
nality relation between normal modes . With .this, additional step 
accurate solutions for higher modes are readily " obtained. 

The basio steps ‘of ‘tile iteration method and a condensed procedure 
for solving problems will he outlined in the following sections. 

Aatual application of the procedure is illustrated 'by a number of 
examples. These examples are intended to serve as a guide In' 'the 
solutions of praotioal problems and at the Bame time to give an 
indication of the simplicity attainable with -this method of analysis. 
All computations are performed mentally or with the aid of a slide 
rule and comparisons made with the few available exact solution 
show that the iteration method gives for praotioal purposes the 
exaot answer. A theoretical verification, of the analysis presented 
is given in the appendix. 


SYMBOIS ... 

■, • <• ... 

( # * . . 1 ' 

L " length of beam; half ‘span for symmetrical beams, full span 

for unsymmetrioal beams 

E v Young's modulus of elasticity 

Gr modulus of elasticity in shear 

I bending moment of inertia 

Ip mass polar moment of. inertia per unit length of beam about 

- J . . 

/ . 

W . f - 

s ' 

m 


axis or rotation 

torsional stiffness constant 

■ . ■ ’ 

weight of beam per unit length 
acceleration due to gravity 
mass of beam per unit length Cl) 
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S ' 

shear 

■* 

M 

■bending moment 

V 

T 

torque ' ■ r 


X • ' ■ - 

distance between stations along beam 


®a 

amplitude of nth mode (n ■ 1,2,3# » • ») 


€ n\ 
and T) n j . 

used in place of to denote very small amplitudes 

t *,*:*■ 


X 

• station coordinate 



elastio spring constant at Jth station. 



concentrated mass at jth station 


*n 

* 

frequency of nth natural mode (bending or torsional) 
vibration, cycles per eeoond 

* 

®a 

• !• 

* . * ' . v .• 

circular frequency of nth natural mode of (bending or 
torsional) vibration, radians per second ( 0 ^ * 2jtf n ) 

w 

^eg 

: :.-td 

equivalent loading 



natural frequency of a spring-mass oscillator, radians per 
second . • • . ; 

* 1 


5 J 

deflection at station , J . 

• 

y or Y 

general notation for deflection 


Y n (0) or 
Yn <0) (x)^ 

i assumed or reasonable approximate deflection of nth mode, 
usually written in terns of a unit tip deflection; 

Y n ^ denotes deflection of a given point; Y n ^^(x) 



denotes modal-deflection function 

\ 

y n or \ 

y n (x) J 

exact deflection of nth mode written in teims of a unit 

tip deflection; y n denotes deflection of a given point; 
y n (x) denotes modal— defleotion function 

- 
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defleotion of nth nod© after i' iterations for both 
derived values of defleotion and values written in 
terms of unit tip deflection; Y^ 1 ' denotes defleotion 

of a given point; Y n ( i )(z) denotes modal-deflection 
function 


l n 


•<i) 


or 


*n 


<*> 


defleotion of nth mode after i iterations where deflection 
is given relative to an axis through center of beam and 
must he corrected to satisfy boundary conditions for 
both derived values of defleotion and values written in 

terms of a unit tip, defleotion; Y^C 5 -) denotes 
defleotion of a given point; Y n *(^(x) denotes modal- 
deflection function 


Yq^*) a constant correction to be applied to Y-^^fx), after 

1th iteration, when a solution is sought for a symmetrical 
mode of a free— free beam 

K^x a linearly varying correction to be applied to Y n , ^ i ^(x) J 

after the ith iteration, when a solution is sought for 
an antisymmetrical mode of a free— free beam 


or 

Yn^(x) j 


deflection of nth mode after all modes lower than nth mode 
have "been removed after ith iteration for both the 
derived values of defleotion and values written in 
terms of unit tip deflection; y n (i) denotes deflection 

of a given point; y n (i)(x) denotes modal-deflection 
function 


$ or $ general notation for rotation in torsional-vibration 
problems (with y and Y replaced by f and $ 
definition of symbols for rotation is similar to that 
for deflection) 


BASIC STEPS OF ANALYSIS 
Iteration Methods 

Bending .— Examination of the differential equation of equi- 
librium for a beam in free harmonio bending vibration 
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d 2 jjj d. ? yn(x) 
cbc 2 dx 2 


*= y n (x) 



shows that a beam vibrating in one of its natural modes lias an 
inertia loading at any point that is proportional to the product of 
the mass intensity and the defleotion at that point. The left-hand 
side of the equation represents the elastic restoring force of -the 
beam and, in the ordinary bending theory of beams, equals the external 
loading. The right-hand side represents the inertia loading which 
may be considered at any -instant of time to be statically applied. 
Deflection functions y n (x} \ which satisfy both oquation (l) and the 
boundary conditions of the beam are oalled the natural modes of the 
beam. The factor is the natural frequency in radians per eeoand 

of vibration of the nth mode. 


A curve approximating the deflection may be assumed and, as 
indicated by the right-hand side of the equation, valuoB proportional 
to the loading at any point may be computed directly by multiplying 
the assumed deflection by the mass intensity at the point. 

The faotor 0^2 can be ignored since it is a constant and sinoo the 
amplitude of the deflection is purely arbitrary. The values of 

then represent the inertia loading at each point. With this 

computed loading, a new defleotion may be -found by any of the known 
methods of engineering beam theory such as the direot— integration 
method; that is, with the loading on the beam as a starting point, 
suooessive integrations give in turn the shear, the moment, the slope, 
and the deflection. 

In the present paper, the process of finding a now deflection 
frcm an assumed or given approximate deflection curve by four 
successive numerical integrations is considered one iteration. In 
the appendix the new defleotion is shown to converge toward the lowest 
mode component present in the assumed deflection. By successive 
iteration, a natural-mode defleotion curve may be found to an arbitrary 
degree of acouraoy. 

An iteration by two different methods for a beam in bending 
vibration is presented in table 1. The problem used for illustration 
is that of a nonuni foxm cantilever carrying a concentrated mass. The 
s u m ma tion method is a simple numerloal-integration procedure to 
determine the shear, the moment, the slope, and the new deflection. 

The second method presented is a form of numerical Integration whioh 
makes use of the ooncept of equivalent concentrated loads. (See 
table 1. ) 
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Most of the steps of an iteration will- become self-evident on . 
inspection of the tabulated computations and the graphical illustrar- 
tions .in figure 1. 

•For the summation method (table 1), the physical constants. I 
and m for the cantilever beam are listed in columns 2 and 3, and 
the initially assumed values of deflection at each station are liBted 
in column h t The beam is divided into equally spaoed stations with 
the root as station 0 and the tip as station 10. In the numerical 
integrations , the distance X between stations may then be carried 
as a common factor as shown. The loading values which are 

proportional to the inertia loading are given in columns 5 and 6. 

The' shear (column 7) and moment (column 8) are then found by the 
•summation process -indicated in figure 2. Except for the initial 
value, the shear. is found by successive addition of the loading 
ordinates. The initial value of shear is found by application, of 
the following simple equation (frcm' fig. 2)j. 

Area = - - ~(2) 


..Addition starts frcm station 10 because the shear is known to be 
zero there. The moment is found by successive addition of the shear 
values] the addition starts at station .10 where the moment is zero. 
Freni the M/l values * (oolumn 9), the slope (column 10) is found 
in the same manner that was used to obtain tho .shear exoept that 
integration prooeeds frcm station 0 where the slope 1 b zero. 
Successive addition of the slope values, starting at station 0, gives 
finally '.the new deflection (column 11). In order to compare the new 
deflection with the assumed deflection, the new deflection is given 
in the next oolumn in terms of a .unit tip ordinate. Burgess's method 
(reference 2) consist's of adding successive station ordinates but 
.a. clear physical interpretation of his procedure cannot readily be 
made . ‘ 

Since the beam has been divided into equal intervals, mass 
concentrations on the beam will not, in general, be looated at 
station points. The procedure, then, is to distribute a concentrated 
mass proportionally to two adjaoent station points on the assumption 
that the mass is supported by a simple beam between the two points. 
When the station points are spaced reasonably close, this local 
redistribution of the loading produces no significant change in the 
’ derived deflection. 

For curves that cannot be reasonably approximated by straight 
lines between stations without resorting to a large number of stations. 
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smother method of performing integrations numerically has been 
employed. In reference 5, the equivalent-load method is presented; • 
this method makes use ,of the properties of a second-degree curve for 
computing the system of concentrated loads at the stations which 
produce the true moment at eaoh station and for computing the concen- 
trated values of M/EH whioh cause the deflection to be correot at 
each -station. The method is easy to apply and gives results whioh 
are quite accurate even when curvatures are appreciable between 
stations . : 

* * ./' i ♦ 

For the equivalent-load method (table 1), the equivalent— loading 
diagram (column h) is oemputed by using equations (a) and (b) given 
in figure 3. These equations • were t alien from reference 5 and are 
immediately applicable only to curves continuous over at least two 
station intervals. In order to keep. mY^(°) (column 3 ) from having 
a discontinuity at station 3, therefore, the concentrated mass 
(oolumn 5) is handled separately-. • Summation of the equivalent loads 
gives the average shear between each station and summation of the 
shear gives the true-moment diagram. The equivalent values of M/l 
(E is oarrled as a common factor) are- then found in the same maimer 
as the equivalent loads with' the exception of the value at station 3. 
The discontinuity in the shear diagram at station 3 oauseB an abrupt 
change in- slope- of the moment diagram. The equivalent M/l value 
at this station is. a amputad by use of equation (c) in figure which 
is derived by applying equation (a) to the ordinates of the smooth 
curve on either side of the abrupt ohange. The deflection (column 11) 
is found after two summations. 

Either the summation method or the equivalent-load method may be 
used for an iteration. The summation method is simple and quick 
and gives good results for the lower modes and frequencies of vibration. 
The method 'of equivalent loads, on the other hand, is particular ly 
suited for higher-mode deterzninati on . For the higher modes, curvatures 
are appreciable in the functions ' ter be integrated, and the greater 
inherent aocuraoy of this method results in better approximations to 
the deflections. , . ( 

Torsion .— The differential equation in torsional vibration is 

«= tc^ 2 I p ^(x) (3) 

This equation is analogous to equation (l) for -bending vibrations. 

The left-hand side of the equation represents the elastic restoring 
force, and the right-hand side represents the inertia .loading. The 


A QJ #£$ 

dx - dx 
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iteration process in torsion is similar to that. in. tending j that 1 b, 
a curve .of angular displacement is assumed - and on multiplication "by 
the mass polar moment of 1 inertia, a measure of the inertia torque 
loading is obtained.. By direct integration anew angular-displacement 
curve is found. Only wq integrations need he performed in ah 
iteration for torsional problems as compared-. with the four integra- 
tions for bending, problems because the differential equation of 
equilibrium-in. torsion is of the second order, whereas the. equation 
in, bending is of the fourth order. ... .. .. .\-jLJl 

A typical iteration for a cantilever in torsional, vibration is 
given in tabular form in table 2 and presented graphically iri figure b. 
Integration. of the torque -loading ordinates (column 5) by the summation 
method gives a value of torque (column 6) midway between '.each 
station. After, dividing ’each torque ordinate by its corresponding 
value of J, another summation gives the rotation (column 8). The 
summation in each integration begins at the station where the function 
being determined. is knowh to ’be zero , Only the summation method of’ 
integration is used in torsional problems because it is not convenient 
to determine equivalent concentrated: inertia torque loads that will 
give the correct internal torque at specific points on the beam. 

Treatment of Various Boundary Conditions 

Bending . - The boundary . oonditionB encountered in vibration 
problems are the .same as those found in problems in statics j that is, 
the beam may have any combination of free, pinned, elastically 
restrained, or fixed ends. In addition, if the beam and its support 
conditions are symmetrical about a center line, tbe equilibrium 
conditions existing at the center line depend on whether the beam 
is vibrating in symmetrical or anti symmetrical modes. For such beams, 
the - two types of modes are found separately by considering the beam 
to be cut. at the center and by applying the^ proper. boundary conditions 
for each type. 

In, the iteration process, the functions determined by numerical 
integration will be correct only if all the boundary conditions are 
satisfied. In the simplest cases, such as a cantilever or a 
symmetrical simply supported ; beam, where the shear, moment, slope, 
and deflection are each known, to' be fcero at specific points on the 
beam, no difficulty is encountered in performing integrations* Sum- 
mations simply proceed' from the. points where the, functions are known 
to be, zero* In other types of beams, there may be more than one known 
boundary condition on some of the functioned and none at all on others. 
For those -functions without boundary conditions, the point at which 
the function is zero’ is generally not knpwn and a, special treatment 
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must, therefore he made ' in the : Integration process* . Although the 
treatment will vary with each' problem,-. -a general method of approach 
'may he outlined* -The method- is to sketch the probable deflection, 
shear, moment, and slope diagrams in* that order. In doing so, the 
conditions that must he satisfied to mak^ each diagram follow from 
the preceding one -will become evident'. It vill he seen that an 
arbitrary constant of integration is introduced when no boundary 
conditions are known for a ; particular function. The constant is 
then carried along and evaluated later when two boundary conditions 
are known for one of the other functions. In some problems, short- 
cut calculations may be made to evaluate the constant immediately. 

As an illustration Of the general approach, consider the -case 
of a free-free beam symmetrical about its midpoint. In figure 5 the 
Series of diagrams on the right represent ' the actual variation of the 
deflection, loading, shear, and so forth, along the half span of 
the beam for the first symmetrical -bending mode. No boundary condi- 
tions are known for the deflection) the shear is zero at both the 
midpoint and. the tip, the moment - - is kero at the tip, and the slope is 
zero at the midpoint. Since an arbitrarily assumed deflection will 
not, in general, produce a loading which will cause the shear to be 
zero at both the midpoint and the tip, the deflection has been 
assumed in two parts as shown in the two diagrams on the left. A 

variable part Y^’ ^ giving deflections relative to the center of 

the beam is assumed, together, with a constant part Y 0 ^ , the 

magnitude of which is to be 'determined later by the shear boundary 
conditions. This separation of the deflection into two parts is 
permissible since the principle of superposition holds for problems 
of this type. • ' ; 

Loadings are computed for both deflections (each deflection is 
multiplied by the mass variation given in column 2 of table 3) and the 
shear diagram for each loading is found by integration from the tip 
inward. In-order to make the shear zero at the center line, the 
ordinates of the S Q diagram must be adjusted to make the center-line 

ordinates of the S' and S Q diagrams . equal in magnitude and 

opposite in sign) the two diagrams then added give the true 
shear diagram S (fig. 5) • If desired, Itr is now possible at this 
point to correct the assumed deflection and to obtain the correct 

loading. Since the value of shear s Q is proportional to Y 0 ^, 

the adjusted’ deflection |U Y ^ added to Y,'^ gives the 

/ f .\ ' . .0 » * * 

deflection Y^ •' When the true shear diagram is known > the moment 
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and slop© diagrams are found without any difficulty. For the new 
deflection integration proceeds frcm the center line which, gives 
the deflections Y^ 1 (1) relative to the o enter of the beam. The 

correction Y 0 (2-) must now be determined. This correction is found 

from Y^/C^.) in the same manner that Y Q ^ was determined from 

Y^*^. After Y 0 W is determined, Y^*^ is corrected to give 
the better modal approximation Yt.( 1). Note that the process of 
determining the corrections Y 0 (-I consists actually of performing 

a large part of another iteration. The method Just described for 
integrating and satisfying boundary conditions is analogous to 
direct analytical integration of the differential equation, where 
the unhncwn constants of integration are carried along until enough 
boundary values are known to permit them to be evaluated. 

A short-cut calculation is derived for the preceding oase from 
the following consideration.' It is recognized that in order to 
prevent translation of the beam the negative and positive loading 
areas must be equal. (See loading in fig. 5«) Expressed mathe- 
matically the condition is 

L 

mY n dx.= 6 (4) 

-or if Y n is broken into two parts, as in figure 5, equation 4 
becomes 



slid- Y q is found to be . 

I nYn’ dx 

- J ° — ( 6 ) 

m dx - 


In equations 4 to 6, and whenever convenient in the equations- to follow, 
the superscripts indicating the iteration number have been purposely 
omitted for simplicity in presentation. The equations are of general 
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form and apply regardless of the number of the iteration. The 
numerator -and .denominator of equation ( 6 ) can conveniently be evaluated 
in terms of , the. station ordinates by means of; the following equation 
which is- derived in the- appendix ;■ 

. - .. * \i ■ ' l 

»* • • . '• * 

A «= ( 0 . 38 a. -I- l»50b q + d + 1.50e + 0.76f + l.^Og 

+ h ? + i +< 1 . 50 j + 0.35k) (7) 

' ‘ , 1 ’■ “ ' 

Equation. ( 7 ) for computing- the "total area under a curve is used 
hereinafter booause of. its simplicity and aocuraoy. This equation 
is based on’ the properties' pf a fifth-degree curve and Is therefore 
applicable to finite intervals in sets of five. After V has been ' 

assumed," Y Q " can be computed directly by Use of equation ( 6 ) and 
then the correct deflection Y n *= Y n ’ — Y Q can be found. This short- 
cut procedure is illustrated by the numerical example in table 3 - A 
complete iteration is given in the section "EXAMPIES". 

In table 4, the general approach is illustrated for the first 
antisymmetric al mode of asymmetrical free-free beam. The diagrammatic- 
presentation of the data is shown in figure 6 . As illustrated in 
the column of diagrams on the right-hand- side the deflection and loading 
are known to be zero at the center of the beam and the shear is zero 
at the tip. The moment must be zero at both the center and tip, and 
no boundary values are known for the slope. The difficulty in this 
set of boundary conditions is that an arbitrarily assumed deflection 
will not produce a loading which will cause the moment to be zero 
at both the oenter and the tip. The moment boundary conditions can 
be satisfied, however, if the assumed ourve is given a proper rotation 
about the midpoint of the beam, A linearly varying correction 

K^^x is therefore integrated along with the. assumed deflection 

Yi’(°), as shown in the two left-hand columns in figure 6 , until 
the moment diagrams are obtained. The two diagrams are added 
(fig, 6 ) to give the true moment diagram. When the correction to be 
applied is known, the corrected assumed defleotion, the correct 
loading, and the shear can be obtained. 

In the next step the -interior point at which the slope is zero 
is unknown, so a zero point is assumed and integration for the slope 
proceeds from this point in both directions. An arbitrary constant 
correction is also assumed. Integration of the incorreot slope 
diagram and the constant correction gives the incorrect defleotion 

■Yj/Cl) and the linearly varying correction K^-^x. The correction 
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K^x to the deflection Y^ r i^ ke' determined’ in the same mariner 

that was used to find the correction E^x to the deflection Y^ 0 ^ 

that is, with Yi»^ and E^x, integration proceeds in the next 
iteration (not shown in the table) .until the' moment diagrams are" 
found, from, which the true correction can he found by the condition 
that the center-line moment must be zero. 

Again an equation can be derived for correcting the' assumed 
deflection before an iteration and the number of integrations is 
thereby reduced. In order, to prevent rotation of the beam about an 
axis at midspan perpendicular to the plane of vibration, the moment of 
the inertia loading about this axis must.be zero. Written. mathemati- 
cally in terms of the two assumed leadings, the condition is 


• .• rl* ... ■ 

/ mCYn 1 n- Kx)x dx = 0 j(8} 

. J o. ' ; 

where K is the slope of a linear c effect ion. Solution for E. 
giveB - . 


(9) 


This equation can be conveniently evaluated numerically by use of 
equation ( 7 ). The use of equation (9) is illustrated by the example 
in table- A complete iteration is given in a subsequent example. 

As an Indication of the flexibility of the general method of . 
approach for • handling boundary conditions, consider next a beam fixed 
at one ’end and simply supported at the other end. The logical 
procedure for performing an iteration on such a problem is given in 
figure 7 / The process should become clear from a study of the sketches- 
and the step-by-step outline of the iteration procedure. 

The foregoing illustrations indicate that rather complicated , 
and even statically indeterminate structures can be handled if they 
are broketi up into their basic component parts. In the examples that 
are presented subsequently; 'the method of adding elastic restraints, 
such as springs, to the beams is shown. The solutions given for these 
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problems majce use of ghpt^-gU^ g<^h|tig^g, similar tp equation (6), 
wMoh affe derived fr© |£sig ghj£che§ suO^af |%£ 
figures 5 §fi§ ■§; ■ * ■; 1( ■ 1 _ - V.- . 

Tcrfiion. - As was trUe iri bbhAingMfibration problem^ no 
difficulties arise' in $h§ fchjjdlirjg gf the ordinary boundary conditions 
for problems in torsion with tnfe exception of the .free— free and 
certain fixed-end beams. For these cases a procedure analogous to 
that in bending must be '.used, 

. * ’ ■ l, • ■ 

* i* , r . . 

• Removal of L'ower-Mode Components • 

Bending .— In the deteriai nation of higher modes by the iteration 
process any components of modes lower than the one to be determined 
must be removed from the assumed ar given deflection. Unless this 
procedure ia followed before each iteration, convergence will tend 
toward the lowest-mode component present ih the assumed deflection. 
Assuming that the shapes of the lower modes are known, the components 
of each mode present in the assumed higher-mode shape can be found by 
use of the orthogonality relation between the normal modes of vibra- 
tion. The equation for computing the amplitude a n of any mode y n 
present in a given deflection curve Y m is 

• r L 

/ lx 

J 0 

Q n - — £ ( 10 ) 

/ lx 

J 0 

This equation is given in reference 1 and, for completeness, a 
derivation is included in the present appendix. ' 

The procedure for removing lower-mode components is illustrated 
in table 5 where the fundamenta 1-mode component is subtracted from 
an assumed second-mode shape fop a cantilever. A graphical repre- 
sentation of thiB procedure is given in figure 8. The integrals in 
equation (10) are evaluated in columns 5 and 8, and the fundamental- 
mode amplitude is computed at. the bottom of the table . The assumed 
second-mode shape minus the component of the fundamental mode is given 
in , column 11 in terms of a unit tip deflection. The same procedure 
applies when n-1 lower modes are being removed from an assumption 

for the deflected shape of the nth, mode. , ‘ 

* 

Torsion .— The method for removing lower-mode components in 
torsional problems is the same as that for bending; 



NACA TN No. 1522 15 


Frequency Determination 

Bending.— Onoe a close approximation for the deflea ti on of a 
given mode has been established, the frequency of vibration can be 
determined. For bending, the frequency may be found from the following 
equation 

V (i) 

(n) 

y n (l + 1) 

where 

y n (i) value of deflection of a point on beam before iteration 

y n ^ + ^ value of deflection of point after iteration when . . 

fundamental mode is being determined; whereas, it is 
deflection of point after iteration and removal 
of lower-mode components for higher-mode determination 

(The deflection y n (i) is usually written in terms of a unit tip 
deflection and, when equation (ll) is used, y n (* + ^ is the 
absolute value of the deflection found in an Iteration. Before the 
next iteration is begun, however, y n ^ + ■**) is for convenience 

written in terms of a unit tip deflection.) The use of equation (11) 
in the frequency determination of a beam is illustrated in a subsequent 
seotion "EXAMPLES." .. .. .. 

The frequency may also be determined from the following 'equations 
which are derived in the appendix; 




^ ay^y^ 1 * ^ cbc 

j 

fmj^ 1 + - 1 -- 1 ] 2 dx 


( 12 ) 


( 13 ) 
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The moment M in equation (12) is the moment which results from the 
loading my n ( i ). Both equations are derived from the energy relations 

of a "beam in. vibration and have been used by other investigators. 

Since "the frequency, as given by these equations, is evaluated from 
the entire defleotion curve, values of deflection whioh "are as accurate 
as those required by equation (11), where the frequency is evaluated 
from the deflections at a point, are not necessary. Thus, if only 
the frequency is desired in a given example, a good value may be 
found from either of these equations by use of a deflection which is 
found from one to two less iterations than those that would normally 
be required if both mode and frequency are to be determined. 


Torsion. — Frequency determination in torsion is analogous to 
that in bending. The frequency may be found from the following 
equation 



where r • 

0^(1) value of rotation of a station on the beam before iteration 

^ value of rotation of station after iteration for fundamental- 
mode determination; value of rotation of station after 
iteration and removal of lower-inode components for 
’• . ‘•higher-mode determination. ’ -- 


SUMMARY OF PROCEDURE 


Thus far, the more important steps in the analysis have been 
explained in some detail. In order-' to facilitate -the actual working 
of problems, the various steps will be summarized in their logical 
order. The procedure will be given for computing the modes and 
frequencies of a beam in bending vibration. An analogous procedure 
is used for torsional problems. 

Fundamental Mode .— The stepB for computing the fundamental mode 
are as follows: 

(l) Divide the beam into a convenient number of equal stations. 
If only the fundamental mode is to be computed, six to eight stations 
will be sufficient. If higher modes are to be determined later, ten 
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stations should he ijsed because any errors present in the fundamental 
mode will have an effect on the accuraoy of the higher modes. With 
the "beam divided into ten stations, a very accurate fundamental mode 
can he obtained. From mass and moment— of— inertia distribution plots, 
read off the values of these variables at each station. Any large 
concentrated masses should be distributed to the nearest station 
points. 

( 2 .) Assume a reasonable approximation to the fundamental mode. 

A convenient way is to sketch the deflected shape of the beam and 
read off the deflection ordinates at each station. 

( 3 ) Perform an iteration with proper regard to boundary conditions 

(4) Use the deflected shape from the previous iteration and 
successively repeat step 3 until the desired degree of convergence is 
obtained. 

( 5 ) Compute the frequency by equation (11). 

Higher Modes.— The steps for computing higher modes are as 
follows: ■- 

(l) Divide the beam into equal stations. For convenience in 
applying: the integration formula (equation (t))» whiah is based on 
a fifth— degree curve, the number of stations should be a multiple of 
five. In order to retain the same degree of aocuraoy for each higher 
mode," the number of stations must be increased when successively 
higher modes are being determined. For the first four modes, however, 
ten stations are sufficient for engineering aocuraoy if .the iterations 
are performed by use of the method of equivalent loads. Determine 
the mass and mament-of— inertia values at each station. 

• (2) Sketch the most probable shape of the higher mode to be 
determined and read off the values of deflection at each station, 

(Each successive higher mode has one additional nodal point.) 

( 3 ) Remove all lower-mode components from the assumed higher- 
mode defleation. 

(4) Perform an iteration with proper regard - to boundary conditions 
and remove loweri-mode components from the derived deflection. 

( 5 ) Successively repeat Btep 4 until the desired degree of 
Convergence is obtained. It may be found that all the lower-mode 
components need not be removed after each iteration. (A more complete 
explanation as to when lower modes ought to be removed is presented 

in the example on the determination of the third mode of a cantilever. ) 
Before and after the final iteration all the lover modes should be 
removed. 
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( 6 ) Compute frequency by equation (11). 


EXAMPLES 

The solution to a number of typical problems is illustrated by 
examples. These examples are presented to serve as a guide when 
actual problems of the type. they represent are being solved and each 
one need, not be studied to understand the basic procedure. given in 
this paper. The examples are: 

Example 1 —Nonuniform cantilever, first bending mode 

Example 2 — Nonuniform aantilever, first torsional mode 

Example 3— Uniform cantilever, second bending mode 

Example 4 -Uniform cantilever, third bending mode 

Example 5 — Free—: free beam with concentrated masses, first 
symmetrical bending mode 

Example 6— Eree- free beam coupled to masses through springs, first 
symmetrical bending mode 

Example 7 — Free-free beam with concentrated masses, first anti—, 
symmetrical bending mode 

Example 8— Beam with concentrated mass and mounted on spring, first 
bending mode 

Example 9 —Beam with concentrated mass and mounted on spring, second 
bending mode 

Examples 1 and 2 are simple illustrations of the iteration 
process for bending and for torsion. Examples 3 and 4 show the 
removal of lower-mode components in the determination of higher modes. 

A comparison is given with the exact solution for modes and frequencies 
Examples 5 to 7 Illustrate the manner in which boundary conditions 
are satisfied for the vibration of free— free beams. Boundary condi- 
tions for beams coupled to springs are illustrated in examples 8 and 9 
and the results are compared with exact solutions. In every example 
one or two complete iterations are shown in detail, and results of 
any other iterations neoesaary for convergence are indicated. 
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Example It Nonuni form cantilever, first bending mod e.- The 
physical constants and first iteration for example 1 are given in 
table 1 and the solution is completed in table 6 . Eaoh step follows 
very closely the procedure outlined for fundnpiental-mode determination. 
It should, be noted that convergence to the fundamental mode in 
cantilever beams is yeiy rapid. After only two iterations the ratio 
was taken at eaoh station of two successive deflections, oolumn 12 , 
and it is seen to be fairly constant. A comparison is given in 
figure 9 between the. assumed defleotion and the calculated funda- 
mental mode. ' 

Example 2; Nonuni form cantilever, first torsional mode .— An 
example of beam torsional vibration is presented in table 7- The 
first iteration is the one used to illustrate the iteration process 
for torsion (table 2). It -is seen in columns 2/7, 9* 3-1* and 14 in 
table 7 that convergence to the fundamental torsional mode for a 
cantilever is slower than for the corresponding mode in bending. The 
rate of convergence in the iteration method is a function of the 
separation of the frequencies of the natural modes and, for a canti- 
lever, the frequencies of the torsional modes are not as widely 
separated as the frequencies of the bending modes. The greater 
number of iterations needed in this particular example can also be 
attributed in part to the large difference in the shapes of the 
assumed mode and derived mode. (See fig. 10.) The labor involved, 
however, in finding a torsional mode is comparable to that in bending 
since only two Integrations are needed for an iteration. 

Example 3? Uniform cantilever, second bending mode A unifonn— 
beam case Is presented in table 8 so tliat a comparison can be made 
between the exact solution and the results obtained by the Iteration 
method. The procedure used is that outlined for higher modes. 
Previously computed ordinates for the fundamental mode are listed in 
column 2 and the assumed second-mode ordinates are, given in oolumn 5 . 
Because the beam is uniform, the solution is somewhat simplified as 
the mass m and moment of inertia I can be ignored during the 
Iterations and are taken into account only in the frequency computa- 
tion. By comparison of columns 35 and 38 it can be seen that the mode 
obtained after two iterations is good enough for most practical 
purposes. If modes higher than the second are to be calculated, how- 
ever, more accurate values should be obtained for the deflection of 
the second mode than those given in oolumn 35 . The results of the 
third iteration are shown in the sketch of figure 11 (a), and they 
present a dose check on the mode and frequency obtained by an exact 
analysis. The denominator of equation (10) is evaluated in column 40 
for use in the third-mode determination, presented in the next 
section. The fact that column 42 sums to zero means that the derived 
second mode 1 b orthogonal to the fundamental mode. 
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Example 4; Uniform cantilever, third bending mod e.- The third 
mode and frequently Deputations for the same uniform team are 
presented in table' 9. .The assumed third mode is given in column 2 
and the fundamental mode and derived second mode from table 8 are 
used in computing columns 3 and p. After each iteration a large 
fundamental-mode component must be removed from the derived deflec- 
tion. It is difficult to remove completely all traces of lower 
modes, and during eaoh iteration the residual lower-mode c exponents 
are amplified; in this oase in the ratios for the first 

mode and (o^/nug ) 2 for the second mode. For a oantilever the 

ratio- (a>g/<«i ) 2 is very large and sifter one Iteration the first- 

mode Impurity has been greatly amplified (column 23). The actual 
amount of amplification of the second-mode impurity (column 2k) is 
small enough to be negligible. Seme labor could have been saved in 
this example by postponing the removal of the second-mode component 
until after the seoond iteration.’ After three iterations and removal 
of the lower modes, the derived third mode (column h 6 ) and frequency 
as oemputed by equation (11) give an indication of the high accuracy 
attainable by the procedure described herein when the beam is divided 
into ten Stations. Figure 11(b) shows the assumed third mode, the 
computed mode, the exact mode, and a comparison of the exact and 
oampiited frequency. 

Example 5 i Free— free beam with concentrated masses, first 
symmetrical bending mode .— In table 10, the first symmetrical bending 
mode and frequency of a uniform free— free beam' carrying concentrated 
masses are computed. A convenient assumption for the deflection is 
one with a zero center-line ordinate (column 5) • In order to insure 
that the shear be zero at the center line, the assumed deflection is 
then oorreoted according to equation 6, the use of which has been 
illustrated in table 3, With the computation starting at the oorreoted 
deflection, given in tonne of a unit tip deflection and multiplied 
by the constant mass of unity (oolumn 8), iteration is performed in 
a straightforward m a n ner. The new deflection (column Ip) is again 
given with a zero center-line ordinate and corrected (column 19) by 
use of equation' 6. Two iterations are sufficient to determine the 
mode and frequency. The computed deflected shape is shown in figure 12. 

Example 6 : Free— free beam coupled to masses through springs j 

first symmetrical bending mode .— Example 6 (table ll) is the same as 
example p except that the concentrated masses are elastic ally o.onneoted 
to the beam by springs. The equation used to make the assumed loading 
(column 5) satisfy the shear boundary condi tions_ls given in figure 13 
and is developed in the appendix. This equation contains terms 
involving the unknown coupled frequency, whioh necessitates an assump- 
tion for this frequency before eaoh iteration. Unless the natural 
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frequency p^ 2 of the spring-maBB system nearly coincides with the 
unlmown frequency a> 2 of the ooupled system, the equation for Y 0 

is relatively insensitive to incorrect assumptions for oj 2. in the 
general case, a more accurate assumption is made for o> 2 before 
each succeeding iteration on the basis of a frequency computed from 
the preceding iteration. 

The system Under consideration in this example involves a 
relatively stiff restraint between the masses and the beam. The 
unlmown frequency Will be only slightly less than the frequency 
obtained in example 5 whioh considers the mass rigidly connected to 
the "beam. In that example (see table 10) the square of the frequency 
was found to be .11.6 "(radians /sea ) 2 ; consequently, an assumption of 
11- ( radians /sea ) 2 was made for thi« case, The initial assumption for 
the deflection is the deflected shape obtained in the last iteration 
in table 10. Once the correct inertia loading is computed (table 11, 
oolumn .8) -iteration -proceeds as' before. Since the assumed deflected 
shape and frequency were very close to their aotual values, one 
iteration results in a uniforia ratio at each station (oolumn 20) 
between the deflection and the assumed deflected shape. The computed 
defleotion is shown as the dashed ourve in figure 12, Elastic 
coupling between the beam and masses is seen to reduoe the frequency 
of the entire system and to cut down the deflection of the oenter 
of the beam relative' to the ends. In table 11 a computation 1 b also 
presented for the force in the spring for the given beam-tip 
deflection. • . . _ 

■Example 7: Free— free beam with concentrated masses, first 

ant lsymaetrioal bending mode .— Table 12 Illustrates the solution for 
• £tie. ’antisymmetrlcal 'modes of the same free— free beam analyzed in 
example 5* A defleotion is assumed (column 7) and by means of 
equation ( 9 ) this deflection is corrected to satisfy the moment 
'boundary conditions. This procedure has been illustrated in table 4. 
With the correct loading (column 12) the integrations are straight- 
forward until the slope (column 19) is to be determined. The point 
of zero slope is unknown and is assumed to be between stations 3 and 4 . 
Integration of the assumed slope diagram results in incorrect deflec- 
tions (qolumn 20).- A true defleotion is computed as, before, however, 
which satisfies the moment boundary conditions ‘and automatically has 
the correct slope at each point. The frequency was aomputed after 
the second iteration. The computed def looted shape is shown in 
figure 14 . _ 

Example 8: Beam mounted on ^pring, first bending mode .- A first 

mode, solution for the problem of a beam vibrating on a spring is 
given in table 13 . The boundary conditions are the same as for the 
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symmetrical free— free 1)6621 modes except that the unbalanced inertia 
loading on the beam must equal the spring reaction. This condition 
is satisfied when the distance Y 0 in figure 15 is computed by the 
equation given in the figure whioh is developed in the appendix. In 
order to solve this, equation an estimation of the probable value of 

CD 2 is first necessary. The following considerations will aid in 
making this estimate. If the spring is infinitely stiff, the lowest 
mode and frequency is the fundamental for the half beam vibrating as 
a cantilever. If the spring has a finite stiffness, however, one 
natural frequency of the system will be lower than the cantilever 
frequency. The more flexible the spring is relative to the beam 
stiffness, the lower the frequency ,of the first bending mode will be 
relative to the cantilever frequency. ' The limiting case is that in 
which the beam is considered to have infinite stiffness and corres- 
ponds to an inelastic mass vibrating on a spring. Call this frequency 

of vibration pj 2 which is given as pj 2 - 2 j 1, The lowest natural 
bending frequency co of 1 the system will always be lower than pj. 

The example worked in table 13 is for the case of a flexible 
beam carrying a heavy concentrated mass and mounted on a relatively 
stiff spring. The oantilever and spring-mass frequencies are first 
calculated in order to form a basis for estimating the probable 
natural frequency. Since pj 2 is several times larger than a> l 2 oarvfc -> 

a natural frequency only slightly lower than the oantilever frequency 

was assumed. After is calculated by use of the formula given, 

the succeeding steps of an iteration are the same as those for the 
symmetrical modes of a free— free beam; The trial frequencies calcu- 
lated in column 17 show that the estimated frequency of 5.9 used in 

calculating Y^^ is as aoo urate an assumption aB could have been 
made. After the second Iteration a few trial calculations for tx^ 
(column 28) showed that an estimated .frequency, of 5.93 gave the proper 
values of Y1 / 2 ). A third iteration is indicated with the resulting 
frequency and mode given in columns 31 and 32.' ■ In figure 16 the 
.deflection is presented with a comparison between the exact and 
computed frequencies. 

Example 9: Beam. mounted on spring, second bending mode .— A 

second-mode solution for the same beam-mass— spring system is presented 
in table 14. The procedure used is that for higher-mode determination. 
The most reasonable assumption for the deflection (column 5) is a 
curve having a unit deflection at the tip and center. The first-mod© 
component present in the assumed ourve is calculated and removed and 
one iteration performed. The deflection (column 16) is given in terns 
of a zero center-line ordinate. This ordinate must be computed by 
use of the equation which was used in the first-mode determination. 
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In order to solve this equation'/ the' frequency must first “be estimated. 
A reasonably acceptable value may he obtained by equating the total 


inertia force ^tu 2 J myg^ d x.j of the assumed loading to the spring 

force (center deflection times spring' constant) and solving, for to 2 . 
This calculation is shown in table 14. The acouracy of to 2 depends, 
of course, on the closeness of. the assumed ourve to the actual 
■"deflection. By use of this va,lue of oi 2 , the equation for Y 0 
figure 15 is solved for Y 0 (l). The value of .Y 0 (l) is subtracted 

from the value of Y 2 1 ^ (aolumn 19) and the first-mode component 
• is' again removed and another iteration perfomed. In eaoh succeeding 
iteration a better value of to 2 iB found by oonsidering the center- 
line shear boundary condition and finally to 2 converges to the fre- 
quency as computed by equation (11). In this example, three ^Iterations 
are sufficient for convergence. The determined mode and frequency 
are presented in figure lo. A comparison is given between the exaot 
and computed frequency. 


CONCLUDING- REMARKS 


The examples presented have shown the manner in which the 
_ determination of vibration modes and frequencies is related to beam- 
deflection theory. By successive approximation, solutions are 
readily found to an accuracy consistent with that to whioh the physioal 
constants in built-up beams are ordinarily known. Although exaot 
solutions, for practical purposes, are obtained .to the differential 
equations, the equations themselves are based on limiting assumptions. 
Since the equations do not include the effects of structural da mp ing, 
rotary inertia, and deflection due to shear, whioh are known to cause 
a ohange in the Bhapes and frequencies of the higher modes, engineering 
Judgement must be used to Interpret solutions obtained from these 
equations . ‘ ' * — " " 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, "Va., July 23, 1947 
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P&termlna.tiQn of the Deflection by the Iteration Method 


Differential equation of equilibrium .- The differential equation 
of equilibrium 6f a beam in free' ha;monia bending vibration is. 

— (m^) = 05^ (Al) 

•. . .d*2\ dcV *: . • 

Tji6‘ left. -Hand 'side of. the equation represents the elastia restoring 
’ forces; of the beam and in the . ordinary bending theory of beams equals 
' the; /eitei-hal. loading. For beams in vibration the loading is composed 
'Of Inertia forces 'which at any instant of time may be considered as 
static forces given by the expression on the right- liand side of the 
equation. Equation (Al) has solutions y « y n and associated 
characteristic values co « c.^, where y n dosoribes the deflection 
of the nth mode of vibration and c% is the frequency of vibration 
of that mode. Except for some simple cases these solutions cannot be 
obtained by, exact, analysis. . . It is -possible* however, to solve 
squat ion,. (-Al) by an. iteration, process in which, successive approxima- 
tions 'found, each approximation more closely representing the 
$rue sqlu^loh.'. ' . . • - - 

■ Integration of differential equation .— Suppose that setae ■, 
.arbitrary, function Y is substituted for . y on the right-hand side 
o^eqi^atiph (Al) ., Except for the factor ofi, which may , be ignored 
hecauao’. y . is purely arbitrary, the equation would be 



This equation can be solved by direct integration. Assume for the 
moment that the solution is, y *= F(x) . From the static loading concept 
indicated in the previous seotion it .follows that a beam with a 

loading oi ri Q iay n would have a deflection y n . A beam loaded with 
forces my n would therefore -have a deflection yn/“n ? '. Thus, if Y 

were an exact solution — say y n — of equation (Al), the solution of 
equation (A2) would be simply 
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y » T(x) = ^ 

<%i 2 


(A3) 


In general, Y is not the exact solution y n , and the proportion- 
ality given "by equation (A3) would not be, constant along the beam. 
Because the process of finding a hew curve from a given curve is a 
converging process, however,, curves givihg a constant proportionality 
may be obtained even though the starting curve is approximate. Hie 
newly found curve is used to determine the loading, and another 
deflection is computed. This .operation is repeated until two suc- 
cessively determined deflection ourves are of constant ratio to each 
other. The final curve found is the fundamental or first mode (n * 1) 
of vibration of the system. This process of finding the curve is 
commonly known as the Stodola method. . The faot. that the process is 
converging and converges to the lowest mode will be shown in the next 
seotion. In the section "Frequency Determination" the way in which 
the method is made to converge to a given higher mode of vibration 
is shown, • ' ... 

Proof of convergence of the iteration process .— The arbitrarily 
ohosen curve Y may be- 'expressed In terms of a series involving the 
normal mode's of vibration of the’ system, timB " 

"• ^°)(x) » a;jyi(x) ’ + a2y 2 (x) + a 3 y 3 (x) + (A4) 


If this series is- substituted for Y on the right-hand side of 
equation (A2) and integration is performed, then, by the principle 
of superposition, the solution will be 

i - ..... 

* (1) (X) - ^Wi(x) + -^%r 2 (x) + -^y 3 (x) + * . , (A 5 ) 


If the prooess is repeated successively with the newly found curve 
as the approximation to the deflection, the curve found after i 
iterations will be 


r(D 


(x) 


a Vi (l) *Te?tf x) * 


“i 1 


2i 


<U2' 


®3 


k y 3 {x) + * 


(A6) 


This equation may be written in the following form 
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r(D 


<X) 


CD' 


2i 


aiyi(x) a^^x). +^^ 2i a 3 y 3 (x) 


(A 7 ) 


Since cd^ < cug, < ^3, . . ., increasing i causes tiie Mgher-jn do 
components to decrease, and the first mode y^ appears more and 
more pure,' The greater the separation of oo^, csg# 03, . . . , the 
stronger will he the convergence to y^. 


Removal of Lower-Mode Components 


The iterative prooess will always cause convergence to the lowest 
mode present in the originally' assumed deflea tion. In equation (A7) 


/ <»1 \ 21 /^l\ 2i / tDj X s1 

Y^n + 1/ < x 0 ^/ < \^n - 1/ 


Not onlyj therefore, is the 


component of each higher mode reduoed during an Iteration but also 
the amount of reduction becomes greater with each successive higher 
mode. Thus, 1 if the components of all the modes lower than the nth 
mode were to be removed completely from a given deflection ourve, 
the solution would neoessarily have to converge to the nth mode 
(the lowest mode remaining). Provided that the Bhape of a given mode 
is known, the component of that mode. in a given curve may be found 
conveniently' by use of the orthogonality of the nomal modes of 
vibration. Suppose that the given curve is expressed by equation (Ah ) 
and that y n is the mode for'Whioh the component is being determined. 
Multiplication through by my n (x) gives 1 


, my n (x)y(°)(x) m a 1 my n ^x)y 1 (x) + asray n (x)y 2 (x) 

* i , 

+ . . . a r ^ay n £? (x) + , , . 

Integration over the length of the 

my n (x)Y(°)( x ) dX a 
UO 

since the y n -tema are orthogonal 
nality condition 


beam gives the relation 


TL 2 

I (x). dx (A9) 

functions defined by the orthogo— 
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my a (x)y n (x) dx = 0 when m n 


(A10) 


In the integrals to follow, the limits of integrations have teen 
emitted "but should he interpreted to he over the half span for all 
beams symmetrical about a center line and over the full span in 
unsymmetrical beams. Solution of equation (A9) for gives 


a n 


/ my n (x)Y(°)(x). ebe 

T, 


(All) 


2 (x) &x 


With known, the nth mode component a^^x) may be readily 
removed from the given curve by subtraction, thus 


y(°)(x) - a E y n (x) = a3jy x (x) + a^r s (x) + . . . 

+ a n~i y + &JX*! y n+ i( x ) + • • • (A12) 


Treatment of Free— Free Beams 

In the case of symmetrical vibrations of free-free beams the 
first normal mode (pure translation, frequency equal to zero) is not 
an elastic mode but must bo included when an arbitrary curve is 
developed in a series involving the normal modes. Thus, the series 
for an arbitrarily chosen deflection Y'(x) is 


Y’(x) =Y 0 + a^Cx) + a^gCx) + a^fc) + * . . (A13) 


where Y 0 is simply a oonetant. Since it is desirable to have the 
chosen curve in terms of the elastic modes alone, the value of Y 0 
must be determined. In order to prevent translation of the beam, the 
inertia loading in each symmetrical mode of vibration must sum to 
zero over the span. The following relation must therefore be true: 
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my n (x) dx = 0 


(A14) 


Multiplication of equation (A13) through ■by m and Integrating over 
the length of the beam by use of equation (Al4) yields 


mY’(x) dx - | mY 0 dx 


(A15) 


which when solved for Y 0 gives 


J'niSr'(x) dx 

*° “ p ■■ 

m dx 


.. (A16) 


Y Q may now' be subtracted from the assumed deflection, thus 


Y(x) = Y'(x)’-Y 0 = a 1 y 1 (x) + + a 3 Y 3 ( x ) + • • • (A17) 


The resulting deflection is. now given in terms of the normal elastlo 
modes of a free-free beam, satisfies the equilibrium condition 
(equation (A15)),and may be substituted in the right-hand side of 
equation (A2). The expression may then be integrated without 
difficulty. J 

For convenience in the. present solution for symmetrical modes of 
free— free beams, the assumed defleotion Y*(x) is always given in 
terms of a zero center-line ordinate. With this assumption, Y 0 then 
represents the center— line deflection. 

In the case’ of anti symmetrical vibrations of free— free beams, 
the first normal. mode is a pure rotation with frequency equal to zero. 
An Arbitrarily assumed curve, representing the deflection of an 
antlsymmetrloal mode and developed in a series involving the normal 
modes, must .include a term corresponding to defleotion due .to 
rotationj thus, . • 


Y’(x) = Kx + a-jy-^x) + agy g (x) + a^y^x) + . . . 


(Al8) 
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where K is a slope and Xx therefore is a linearly varying 
defleotion. In order to prevent rotation, the moment of the inertia 
loading about the midpoint, of the beam must eq.ua! zero. Thus, 


n 


J xm y n (x) clx - 0 


(A19) 


Multiplication of equation (Al8) through by am and integrating 
over half the length of the beam by use of equation (A19) yields 


* P 

xnff*(x) cbc * j Kx 2 ® cbc 


J 


which, solved for K, gives 


E 


xnff 

y 


xaK'(x) dx 


1 


X%L 


fix 


(A20) 


(A21) 


After this value of K is multiplied by x, Sx may be subtracted 
from the assumed deflection; thus. 


Y(x) * Y’(x) - Ex * a iyi(x) + agygCx) + + * • * (A22) 

This expression may now be substituted in the right-hand side of 
equation (A2) and integrated without difficulty. 


Treatment for Beams Coupled to Springs 

In order to analyze structures which contain an elastic restraint, 
such as a spring, the equilibrium of forces existing between the 
structure and the restraint must first be considered, For the case 
Bhown in figure 15, a deflection cannot be assumed directly because 
the equilibrium condition (spring force equals unbalanced shear at 
center line) may not be satisfied. On assuming the variable part 
Y*(x) of the deflection, however, the constant part Y 0 may be 
determined by equating the total inertia load to the force in the 
spring; thus, 

’’1 

<0 2 nff(x) 6x = -Ctf 0 (A23) 

Vj 



30 


NACA TN No. 1522 


Since Y(x) has been "broken u.p /into, the -variable part Y'(x) and 
the constant Y q -, equation (A23) nay be written 


or 


(*> T o] 


dx 


-ofY, 


(A ?k) 


Solution- '-for Y Q gives • 

J* iriY' (x) dx 

• y q = r : 

■ - + P m dx 

a> J 


(A 25) 


The generalization of equation (A25) for a spring placed at the Jth 
station along the beam is 


Y„ w 
o 


7 


n k'(x) dx 


_ 6 J g ,1 


or 


hi- 


dx 


(A2 6) 


where 8 .« is the deflection of the beam at the jth. station relative 

V 

to the center of the beam and ctj is the elastic constant of the 
spring at the jth station. Subtraction of Y Q from Y'(x) then 
gives 


Y(x) « Y'(x) - Y q * a^y-^x) + agy 2 (x) + a ^(x) (A27) ' 


■ which is the .arbitrarily assumed curve corrected to satisfy the shear 
boundary conditions. " * ; 1 


For the case of a free-free beam" "coupled, to spring-mass 
oscillators (fig. I3) , "the total inertia load must again equal the 
force in the spring . At any instant bf time 


o> 


7 


ifif(x) dx + F = 0 


(A26) 


where F is the sprang force.. For convenience, the equilibrium of 
forces at maximum displacement will be considered. If equation (30) 
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in referenoe 1 (p, 59) is multiplied through by a.ja 0 and then 
ajy 0 is replaced by W, a 0 by (-Y 0 + 8 ^), and by p^, the 
following expression is obtained for the spring force IF: 



(A2 9) 


Since the natural frequency pj of the simple oscillator is given by 



equation (A 29 ) may be written 


(A30) 



(A31) 


Stibstituting this value of F in equation (A28) and writing Y(x) 
in' terms of Y’(x) and Y 0 gives the expression 



(A32) 


(A33) 
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The generalization of equation (A 33 ) for any number of elastically 
mounted masses symmetrically placed at J points along the beam is 


PiriY’Cx) dx + X" — 

J / 1 _/jlV 

4~ W 

Y o 


(A34) 


I 


m dx + 




1 


1 


Subtracting the value Y q from the variable deflection Y’(x) gives 


Y(x) = Y*(x) - Y 0 - a^y^x) + a^gCx) + a 3 y 3 (x ) + 


• • # 


(A35) 


which represents an arbitrarily assumed deflection of the beam 
corrected to satisfy the equilibrium condition between beam and 
spring. 


Frequency Determination 

In equation (A 7 ) all terms except the first become negligible 
when 1 is large enough, and the equation reduces simply to 

Y^(x) = -i — a^Cx) (A36) 

w l 2i 


If one more Iteration had been performed the equation would have 
reduced to 


Y (l+ 1 ) (x) 


^(1+D 


a i 7 l( X ) 


(A37) 


Division of equation (A 36 ) by equation (A37) gives for the square 
of the frequency 
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“l 2 


Y (i) 

Y (l+l> 


(A38) 


Thus, after enough iterations have "been performed to cause the higher- 
mode components to "become small, the ratio- of two successively found 
defleotions is the square of the natural frequency of the fundamehtal 
mode of vibration. _ 

• Although equation (A38) denotes the frequency of the fundamental 
mode, a similar equation can "be derived for the higher modes. Suppose 
that the nth mode is being determined and that i iterations with 
the neoessary removal of lower-mode components have been performed; 
then, with the removal of all modes lower than the nth mode after the 
ith iteration, the resulting ourves may be expressed by the equation 

yn^Cx) = exyqCx) + + • • • 


+ * * * + • • • 


(A 39 ) 


The small coefficients eg* up to remain beoause in the 

numerical procedure each mode lower than the nth mode cannot be 
removed precisely. These values, however, are very small in comparison 
with values of a n . Iterating once again, beginning with equation (A39), 
gives the curve 


Y n ( i+1 )(x) = — yqfe) + *~~*y 2(x) + . . . -f£.y n (x) 


o>l 

+ « * 


0*2 


“n c 


u. 


• — + • • • 

V p 


(A 4 o) 


If 


is factored out of the right-hand side of equation (A 40 ), it 




is seen that, in comparison with the first-mode component has 

been amplified by the factor > the second mode by the factor 

W ) 
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K 2 

\“2, 

nth mod© gives then the deflections 


:) 


and so forth. Removal again of all modes lower than the 


"n( i+1 )( x ) 3 wiM + + • • • y n < x > 

+ . • . — ^ yp(x) • • • (A41) 

Oip 

a~ 

where q, tip are extremely small in comparison with — « . 

. 

Division of equation (A 39 ) by equation (A4l) and neglecting all 
terms in • € > t*, and q gives } as a close approximation to the square 
of the frequency, the relation 


«n 8 


yjjCi+i) 


(A4S) 


Alternate Method of Frequency Determination 


Another method for .determining the frequency can be found from 
the energy expressions for a beam in vibration. It can be shown that 
the potential energy of bending U for a beam in its maximum 
displaced position is given by the expression 


/ g \ 

Us- | 3311 ( 

\dx 2 . 


dx 


(M3) 


and the kinetic energy of the beam V as it passes thorough the 
equilibrium position is 


T rt 


cc^my 2 dx' 

These energies must he equal, hence 

,2 


1 

•o 




dX <a i 


n o 

CD < TTiy c dx 


Solution of this equation for jx£ gives^ 


or 


J* my 2 dx 


dx 


(Akh) 


(A45) 


(A46) 
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With, the use of the relation M 
written 


in- 


equation (A46) may he 




/* 


my 2 dx 


(A47) 


As applied to the iteration process, the moment M in equation (A47) 
is associated with the deflection that is derived. That is, y in 

equation (A47) corresponds to y(i+l), the deflea tion derived in an 
iteration, where y(i) is the given deflection or the defleotion 
used at the start of the iteration. By use of equation (A42), the 

frequency may he written in terms of y(i) and the moment rather 
than y(i+l) and the moment. 

Substitution of the expression for y(l+l) from, equation (A42) 
for y in equation (A47) and solving for co^ 2 gives 



(A48) 


Thus, after a reasonably good approximation to the deflection of a 
mode has been established, the moment which results from a loading 
computed from this defleotion is found, and by use of equation (A48) 
a good value of the frequency of that mode a an be determined. 

Equation (A42). expreBBes the frequency in terms of the 
successively found ordinates at any station along the beam. 

Equation (A46) can be transformed so that the frequency is expressed 
in terms of all the ordinates along the beam. The transformation 
is given herein as a matter of interest. On integration of the 
numerator by parts, equation (A46) would appear. 



(A49) 



3^ 
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For the (i + 1) deflection this equation would he 

dx 


(Si 


Byd+l) jL 

Pr d'VU+Dl 

ax 2 

dx 2 


m 


V {1+1) ’ 


(A50) 


dx 


If y(l+l) were the exact deflection. It would necessarily have to 
satisfy equation (Al) j thus, 




*2 fl+1) 
EE(x) . d JL - . ... 
dx^ 




(A51) 


With the use of equation (A42), equation (A51) may he written 


_df 

dx £ 


El(x) 


d^(^) 


dx £ 


5= m(x)y^ i ^ 


(A52) 


Substitution of equation (A52) in equation (A50) then gives 
for the frequenoy the relation 


cn' 


,2 


JZy( i )y( i+1 ) dx 
J'm [y^ i+1 ^ ] dx 


(A53) 


This equation' gives, in effect, a weighted average of the frequencies 
given by equation (A42) for all the points on the beam. The accuracy 
of the frequency obtained from equation (A53) will in general be 
greater than that obtained from equation (A42). Any local error in 
defleotion will cause no appreciable error in the frequenoy as given 
by equation (A53)j whereas an error in the deflection of a point will 
cause a like error in the frequenoy if determined by equation (A42) 
with the deflections at that point. 


Iteration as Applied to Torsional Problems 

For torBional problems the differential equation of equilibrium 
is 
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_ _d_ 

&x 

This equation is analogous to equation (Al) for tending vibrations . 

The left-hand side of the equation represents the elastic restoring 
forces and the right-hand side represents the inertia loading. A 
procedure similar to that used in the. case of bending vibrations is 
used to solve equation (A5 4 ); that is, a curve representing the mode 
to be determined is assumed and is substituted on the right-hand side 
and, by direct integration, a new curve is found; and so on. Only 
two integrations need be performed in an iteration for torsional 
problems as compared with the four for bending problems because the 
differential equation of equilibrium in torsion is of the second order, 
whereaB the equation in bending is of the fourth order. The proof of 
the prooees for obtaining a solution by Iteration follows closely ■ 
that given for bending vibrations and, therefore, no further 
discussion will be given. 


GJ M 
cbc 


= 


(A54) 


Derivation of Equation Used in Numerical Evaluation of Areas 

Suppose the area under a given curve is to be determined. In 
general, any part of the curve over a finite' interval may be approxi- 
mated by a second—, third-, or higher-degree . curve ; the accuracy of 
approximation, of course, increases with the degree of the curve. 

In the present analysis a fifth-degree curve has been used for 
convenience. 

Consider the plot of a curve y = f (x) . Let z = ~ in order to 

X 

obtain e hLmansionles s coordinate, , and let tne ordinates y at 
z = 0 ; "... <£_- 3. 4, 5 be a, b, o, d, e, f, respectively. It oan be 
readily verified that the following equation represents a factored 
general fifth-degree equation having the required values 
y — a, o, c, . . .' at z — 0 , 1 , 2 , . » *■ 



oo 

CO 


y = - - !)( 2 - 2)(z - 3)(. ~ *0( z - 5) + ^z( z - 2)(z - 3)(z - k)(z - 5) 

- ~c z (z - l)(z - 3)(z - ^)(z - 5) + r~dz<z - l)(z - 2)(z - h)(z - 5) 


_ iez(z - l)(z - 2)(z - 3)(z - 5) + -^fz(z - l)( z - 2)(z - 3)(z - k) 


(A55) 


The area under the curve in the interval z = 0 to z = 5 is found simply by integration of 
this equation; thus. 


P5X 


A = 




y dx = X 


00 


00 


y dz = (0.38a + 1.50b + c + d + 1.500 + 0.38d) (A56) 

lw 


By application of equation (A56) to adjacent sets of five intervals the complete area under 
any curve may be found. Usually ten intervals are adequate for most purposes; the area in 
terms of the eleven coordinates would then be 


A » (0.36a + 1.50b + 0 + d + 1.50e + 0.76f + 1.50g + h + i + 1.50,) + 0.38k) (A57) 


For most purposes an area may be evaluated most conveniently in tabular foim. Columns 3 and 6 
in table 3 are Illustrations of such an evaluation. 
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TABLE 2.— ILLUSTRATIVE EXAMPLE OP ITERATION PROCESS FOR A 
CANTILEVER IN TORSIONAL VIBRATION 

[l « 100 in.} X - 10 in.; 0 =. 4,000,000 psi] 


Station 



tOX=L ► 

■ « ■ « ■ i.i 

3 4- S 6 7 & 9 iO 


1 

2 

3 

4 

5 

6 

7 

8 

9 

Sta- 

tion 

H 

I 

P 

*L<°> 

Vi (0) 

T 

T/J 



Common factors ■■ .... - y , 

X 

X 

\ 2 /d 


10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 

90 

145 

235 

380 

595 

875 

1255 

1750 

2360 

3000 

120 

135 

163 

202 

259 

334 

471 

726 

1100 

1480 

1850 

1.00 

.90 

.80 

.70 

.60 

.50 

.40 

.30 

.20 

.10 

0 

120.0 

121,5 

130.3 

141.3 

155.2 
167.0 

188.2 
218.0 
220.0 
148.0 

0 

* 60.2 

181.7 
312.0 
453.3 

608.5 

775.5 

963.7 

1181.7 

1401.7 

1549.7 

0 .669 
1.252 
1.329' 
1.192 
1.023 
.886 
.768 
.676 
.594 

.516 , 

’ 

8.905 
8.236 
6.984 
5.655 
4.463 
3.440 
2.554 
1.786 
1.110 
^ .,516 
■' 0 

1.000 

.925 

.785 

.635 

.501 

.386 

.286 

.200 

.125 

.058 

0 


a 6o 2 „ (3.x 1201 + 121.5 
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TABLE 3.“ SATISFACTION OF BOUNDARY CONDITION FOR 
THE FIRST SYMMETRICAL MODE OF A FREE-FREE 
BEAM BY A SHORT-CUT CALCULATION 


1 

2 

3 

4 

5 

6 

7 

8 

Sta- 

tion 

m 

Z 

Ti'CO) 

nffl* (°) 

z 

Yi(°) 

(a) 

Y X ( 0 ) 

Common 125* 

factors 244 



125* 

144 



10 

10 

3.8 

1.00 

10.0 

3.8 

O .806 

1.000 

9 

13 

19.5 

.84 

10.9 

16.4 

.646 

.801 

8 

16 

16.0 

.68 

10.9 

10.9 

.486 

.603 

7 

20 

20.0 

.53 

10.6 

10.6 

.336 

.416 

6 

25 

37.5 

.40 

10.0 

15.0 

.206 

.256 

5 

31 

22.6 

.29 

9.0 

6.8 

.096 

.119 

4 

38 

57.0 

.19 

7.2 

10.8 

- .004 

-.005 

3 

47 

47.0 

.11 

5.2 

5.2 

- .084 

-.104 

2 

59 

59.0 

.05 

3.0 

3.0 

- .144 

-.179 

1 

75 

112.5 

.01 

.8 

1.2 

- .184 

-.228 

0 

100 

38.0 

432.9 

0 

0 

0 

83.7 

- .194 

-.240 


a Y = Y *(0) _y (0) 

11 o 


vhere 


Y (0) - 2 nQr l ,(0) 

O rm 


g3--7- 

432.9 


0.194 
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TABLE 4.- SATISFACTION OP BOTJHMHT CONDITIONS FOR THE FIRST A HTIffiM IBT R ICAL 
MODE OF A FREE-FREE BEAM BT A SHORT-CUT CALCULATION 


1 

2 

3 

* 

5 

6 

7 

8 

9 

10 

11 

Sta- 

tical 

a 

X 

ax 2 

£ 


snT 1 "( 0 ) 

£ 

K«»x 

Tl<°) 

(a) 

T^O 1 ) 

Cannon 



r 2 

125XL 2 







factors 

Jj 

I» 

TW" 


L 

144 




10 

10 

1.0 

10.00 

3.80 

l.QO 

10.00 

3J30 

-0.418 

1.4-lB 

1.000 

9 

13 

.9 

10.52 

15.80 

.60 

7.01 

IO .51 

- .376 

.976 

.689 

8 

16 

.8 

10.23 

10 .,23 

.24 

3.07 

3.07 

- .334 

.574 

.405 

7 

20 

.7 

9-86 

9.80 

! ~.10 

-1.40 

-1.40? 

- .293 

.193 

.136 

6 

25 

.6 

9.00 

13.50 

~-38 

-5.70 

-8.55 

- .251 

-.129 

-.091 

5 

31 

.5 

7.75 

5.89 

j ~.62 

- 9 . 6 a 

-7.29 

- .209 

1 

• 

iff 

H* 

-.290 

* 

38 

.4 

6.08 

9.11 

j -.70 

: -10.62 

- 15-92 

- .167 

-.533 

-.37 6 

3 

47 


4.23 

4.23 

-.60 

-8.46 

-8.46 

- .125. 

-.475 

-.335 

9 

59 

.2 

2.36 

2.36 

-.42 

-*.96 

-4.96 

- .084 

-.336 

-.237 

1 

75 

,X 

.75 

1JL2 

- ,22 

-1.65 

'2.^ 

- .o4e 

-j .78 

-.126 

0 

IOC 

0 

0 

0 

0 

0 

0 

0 

0 

0 






75.84 



-31.67 

- 


* 


Ti (0) vhane e(0> 


L E 




75.84 * v*> 


— 0 , 00 * 18 . 
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TABLE 5._ REMOVAL OF THE FUNDAMENTAL-MODE COMPONENT II? AI? ASSUMED SECOND 
MOTE FOR A CANTILEVER IN BENDING VIBRATION 




1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

tion 

m 

7l 

n(yi) 2 

E 

Tg(0) 

myiT 2 (°) 

z 

a 1 (°)y 1 

y 2 (°) 

(a) 

y 2 <0) 

Common factors 

^ 

125 X 

144 



1254 

144 




10 

1 

1.000 

1.000 

0.380 

1.00 

1.000 

0.380 

-0.514 , 

1.514 

1.000 

9 

1 

.847 

.718 

1.077 

.50 

.423 

.634 

-.435 

.935 

.617 

8 

2 

.698 

.974' 

• 9^4 

.06 

.084 

.084 

-.358 

.418 

.276 

7 

3 

.545 

.891 

.§91 

-.32 

-.524 

-.524 

-.280 

-,o4o 

-.026 

6 

4 

.409 

.670 

1.005 

-.54 

-.884 

-1.325 

-.210 

-.330 

-.218 

5 

5 

.287 

.412 

.313 

-.60 

-.861 

-.655 

-.147 

-.453 

-.299 

4 

6 

.180 

.194 

.291 

-.52 

-.561 

-.841 

-.092 

-.428 

-.282 

3 

7 

.105 

.077 

.077 

-.36 

-.264 

-.264 

-.054 

-.306 

-.202 

2 

8 

.044 

.015 

.015 

-.18 

-.063 

-.063 

*’ -.023 

-.157 

-.104 

1 

9 

.011 

.001 

.002 

-.04 

-.004 

-.006 

-.006 

-.034 

-.022 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

■ 




5.025 



-2.580 





a - Yg^ ~ a l^°^ y l vh ® 1 ' 8 


W7l T g^°^ _ -g.580 
Ea(Ti) 2 5.025 


-0,514 
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TABIE 6.- HOHDHIPORK CANTILEVER, FIRST BKHDUW BODE 


ra 


3 

■a 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Sta- 

tion 

Ti(D 

«1 (1 > 

r eq 

P *1 

3 

M 

MA 

(«A)eq 

Slope 

ro 

T x ( DA <a > 


Cooton fan torn- — >• 

XA2 

\A 2 

xAa 

X 2 A2 

X 2 A2 

X 3 A*4 

xVi44b 

xVi44e 



10 

1.000 

10,00- 

60.4 



0 

0 



26221 

5.49 

1.000 






60.4 




3833 




9 

.853 

10.23 

122.9 



60.4 

2.74 

37 


22388 

5.49 

.853 






183.3 




3796 




8 

.706 

10.59 

126.3 



243.7 

9.74 

118 


18592 

5.48 

.709 



" 



309.6 




3678 




7 

.566 

10.19 

123.0 



553.3 

17.83 

214 


1*914 

5.47 

.569 






432.6 




3464 




6 

.43* 

9.55 

113.8 



985.9 

25.95 

311 


11*50 

5.46 

.436 






564.4 




3153 




5 

.31* 

8.16 

97.5 



1532.3 

33.30 

398 


8297 

5.45 

.316 






643.9 




2755 




4 

.210 

6.30 

75.5 



2176^2 

38.85 

466 


55*2 

5.45 

.211 






719.* 




2289 




3 

.123 

*.31 

51.7 

443 

1214.1 

2895.6 

43.90 

2fj| 

1753 

3253 

5.45 

.124 

2 

.0566 

2.26 

27.6 



4109.7 

53.40 

639 


1500 

5.44 

.0572 






12*1.7 




1114 




1 

.01*6 

.66 

8.9 



5351.4 

60.85 

728 


386 

5.44 

.0147 






1250.6 




386 




0 

0 

0 

0 



6602.0 

66.02 

386 


0 

: 

0 


m 2 - 5.45 (radians/ seo) 2 
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TABLE 7 .- 50HDBIF0BM CANTILEVER, FIRST TORSIONAL MODE 


1 

2 

3 

4 

5 

6 

7 

& 

9 

10 

11 

12 

13 

14 

Sta- 

tion 


Vi (1) 

T 

T/J 

0 l ( 2 ) 



V 3) 

h ay 



h (5) 


*<3) 

Conmon factors- ^ 

X 

X 

x 2 /a 


x 2 /o 


x 2 /g 


X 2 /t! 



10 

1.000 

120.0 

60.6 

0.674 

8.314 

1.000 

8.068 

1.000 

7.947 

1.000 

7.919 

5050 

1.000 

9 

.925 

124.9 

185.5 

1.280 

7.640 

.919 

7.395 

.916 

7.276 

.915 

7.248 

5050 

.915 

8 

.785 

128.0 



6.360 

.765 

6.123 

.759 

6.007 

.755 

5.979 

5050 

.755 




313.5 

1.333 










7 

.635 

128.0 

441.5 

1.163 

5.027 

.605 

4.805 

.595 

4.697 

.591 

4.671 

5060 

.590 

6 

.501 

129.8 

571.3 

.962 

3.864 

.465 

3.644 

.451 

3.570 

.449 

3.549 

5060 

.449 

5 

.385 

129.0 

700.3 

.800 

2.902 

.349 

2.715 

.336 

2.654 

.334 

2.638 

5060 

.333 

4 

.286 

134.8 

835.1 

.66 6 

2.102 

.252 

1.950 

.242 

1.903 

.240 

1.891 

5070 

.239 

3 

.200 

145.2 

980.3 

.560 

1.436 

.173 

1.323 

.164 

1.288 

.162 

1.280 

5060 

.161 

2 

1 

.125 

137.3 

1117.6 

.474 

.876 

.105 

.802 

.0995 

.779 

.0980 

.774 

5060 

.0976 

1 

.058 

85.9 

1203.. 5 

,402 

.402 

.048 

.366 

.0454 

.355 

.0446 

.353 

5050 

.0445 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 


0 


m 2 ■ 5050 (radians/ see) 2 
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TABOUB 8.- SECOHD EHJTDIW3 NODE OP A mUFQKM CAjrmUVH? - Continued 
in; A - 10 ln.j m = 5 I = 50 la. 4 } E - 10,000,000 pal] 


17 18 


Station I 2 ^ yjTgW 


20 21 


>V 125X5|i 2 A 11 * A*tt 


144%! 144EX 144EI 


23 24 


■fVj T S (1) T 2 {1) p eq 


12 I 12 


2 6 27 


K eq Slop# 7 2 (2) *1 T 2 (2) 


»l (2 Vl 


■A3 ba 3 nA 4 I A 4 ® 2 I 125A5* 2 A 4 m 


144BI 144®I 144®I 144®I 144 £ ®I 144EI 


1.000 

.519 

.0625 

t323 

-.589 

-.703 

-.665 

-506 

-.286 

-.0872 

0 


5.026 

6.252 

.821 

-3.757 

-6.916 

-8.284 

-7.859 

- 6.011 

-3.453 

-1.158 

0 


5.03 

11.28 

12.10 

8.34 

1.42 

- 6.86 

-14.72 

-20.73 

-24.18 

-25.34 


0 

5.03 

16.31 

28.41 

36.75 

38.18 

31.32 
16.60 
-4.13 
-28.31 
-53.65 


449-9 

368.0 

193.2 

-53.0 

-340.9 

-270.9 


(l) - -22L. - 82.3 

a- (2) „ r-3 3 -?,., _ -in 


2.879 

2.879 

HAH ORAL ADVISOR! 
CMtllWSS FOR AERONAUTICS 
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TABLE 8.- 3BC0HD BKHDINQ MODE OF A TJHIFORM CABTIIEVER - Concluded 
[l - 100 la.; X =■ 10 In.; ■ - 5 1-50 in. 11 . e - 10,000,000 poi] 




32 

33 

34 

35 

36 

37 

38 

39 

40 

4l 

42 



Ts( 1 )/T2^ 2 J 

7 S < 2 > 

7 2 (3) 

72(2)/ 7i (3) 

W 

y a (3) 

(J 2 ) 2 

£ 

Tl7 2 

£ 

Canon 

ml* 



mA* 1 




125 An 


125Xm 

factors 

144EI 



144EI 




144 


144 

10 

2927 

491 

1.000 

2959- 

486 

1.000 

1.0000 

0.3800 

1.000 

/•/ 

0.380 

9 

1531 

488 

.523 

1550 

486 

.523 

.2740 

.4110 

.451 

.677 

8 

201 

448 

.0686 

207 

478 

.070 

.0049 

.0049 

.051 

.051 

7 

-931 

500 

—318 

-938 

488 

-.317 

.1005 

.1005 

-.187 

-.187 

6 

-1726 

491 

-.589 

-1744 

486 

-.589 

.3470 

.5210 

-.272 

-.408 

5 

-2085 

485 

-.713 

-2112 

486 

-.714 

.5090 

.3870 

-.243 

-.185 

4 

-1994 

480 

-.681 

-2022 

485 

-.684 

.4675 

.7010 

-.157 

-.235 

3 

-J.533 

475 

-.524 

-1557 

485 

-.526 

.2765 

.2765 

-.072 

-.072 

2 

-876 

470 

-.299 

-891 

484 

-.301 

.0906 

.0906 

-.019 

-.019 

1 

'-269 

467 

-.0919 

-274 

484 

-.0926 

.0086 

.0129 

-.0016 

-.002 

0 

0 

— 

0 

0 

— 

0 

0 

0 

0 

0 









2.8854 


0 
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CO 


^ • ■■■ » 486 (radians/a®« 2 ) 


72 


(3) 
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a y (°) _ T (0) _ - (0)_ _ - (0) 
y 3 3 *1 T l *2 


vhaM * a ' 0) - T ^'-°- 0785 

md «„(o) - 0.0338 
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TABLE 9.- THIRD BEHDIHD HDDS 09 A UKFORH OABTHJffVER - Continued 

» loo ln.j x - 10 in.; a - 5 *Ai52£l^.j I - 50 la/j 1 - 10,000,000 paij 


17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 30 

31 

32' 

Station 

T 3 (l) 

.ft” 

z 

VjW 

Z 

*1^71 

a S (l) 7 8 

column 23 + 
COlUML 24 

r 3 {1) 

t 3 (i) 

a 

3 

n 


Slops 

Com on , aX^ 

1MBI 

125A 5 n 2 

xV 

144*1 

125x5a 2 

l44*I 

aX* 

nx“ 

mX* 


±n 

12 

A. 

12 

& 

12 

X^B 

144SI 

X^a 

factors 

lW^I 

144 2 EI 

144H 

144ei 

144*1 

144*1 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 

. 

586.6 
302.8 

69.3 

-44.4 

-17.6 

103.0 

227.6 
268.8 
200.3 

72.5 

0 

586.6 

260.5 

50.3 

-26.2 

-8.1 

35.2 

52.2 

36.7 

12.8 
1.2 
0 

223.0 

376.0 

50.3 
-26.2 
-12.2 

26.7 

78.4 

36.7 

12.8 
1.8 
0 

586.6 

158.0 

4.8 

14.1 

10.4 

-73-8 

-155.2 

-141.7 

-60.2 

-6.7 

0 

223.0 

237.O 

4.8 

14.1 

15.6 

-56.0 

-232.8 

-141.7 

-60.2 

-10.0 

0 

267 

230 

194 

158 

123 

91 

61 

36 

17 

4 

0 

-a 

-1 

0 

1 
1 
1 
1 
1 
1 
0 
0 

265 

229 

194 

159 

124 

92 

62 

37 

18 

4 

0 

322 

73 

-125 

-003 

-142 

11 

165 

232 

182 

69 

0 

1.000 

.226 

-388 

-.630 

-.441 

.034 

.512 

.720 

.565 

.214 

0 

4.371 

2.872 

-4.284 

-7.129 

-5.006 

.411 

5.874 

8.277 

6.584 

2.705 

4.371 

7.243 

2.959 

-4.170 

-9.176 

-8.765 

-0.891 

5.386 

11.970 

14.675 

0 

4.37 

11.61 

14.57 

10.40 

1.22 

-7.54 

-10.43 

-5.04 

6.93 

21.60 

55.3 

135.0 

167.7 

119.8 
15.1 

-84.6 

-116.9 

-53.9 

85.9 

98.9 

422.3 

367.0 

232.0 

64.3 

-55.5 

-70.6 

l4.o 

130.9 

184.8 

98.9 

765.3 

-6.2 


*i(x) 


15^2. « 267 
e.879 



- 6.2 

2.B85 


-2.15 


imOlil.ADTTflOai 
C d W NR JfflROWWnoa 
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TABLE 9.- THIRD BEHDHC HOE® 0? A UHIFUHH CAmtHEVBR - ConoludeO 
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TABLE 10 FIRST SYMMETRICAL BEHDIHG MODE OF A FRfeE-FREE BEAM CARHYIHG COHCEHTRATED MASSES - Concluded 
[l = 100 In.; X - 10 In. j 1-5 in. 4 ; E - 10,000,000 pal] 


18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

Station 

^d) 

Ti(°)> x d) 

xr x <« 

F eq 

F eq 

S 

M 

^eq 

Slope 

*l’( 2 > 

Tl ’(2) 

£ 

T x (2 ) 

Tidy^d) 

Tl(2) 

Common 

factors 7 

X 4 



_X_ 

12 

X_ 

12 

_X 

12 

x£ 

12 

X3 

X3 

X 4 

mX 4 

125x5m 

X 4 



144EI 

l44Et 

144EI 

144EI 

144EI 

144?EI 

144ei 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 

59240 

47860 

36560 

25510 

14980 

5280 

-3230 

-10190 

-J.5290 

-18390 

-19430 

12.1 

11.9 
11.5 

10.9 

10.4 
10.0 
12.7 

11.5 
11.3 
11.5 
11.3 

l.ooo 

.808 

.616 

.*31 

.253 

.0892 

-.0545 

-.1720 

-.258 

-.310 

-.328 

5.61 

9.70 

7.40 

5.18 

3.05 

1.09 

-.62 

-2.03 

-3.06 

-3.68 

-1.95 

-20.65 

5.61 
15.31 
22.71 
27.89 
30. 9* 
32.03 
31.41 
8.73 
5.67 
1.99 

0 

5.6 

20.9 

43.6 

71.5 

102.5 

134.5 
165.9 

174.6 

180.3 

182.3 

77 

258 

529 

861 

1231 

1613 

932 

1013 

2092 

2160 

1090 

11856 

11779 

11521 

10992 

10131 

8900 

7287 

5342 

3250 

1090 

82148 

70292 

58513 

46992 

36000 

25869 

16969 

9682 

4340 

1090 

0 

82148 

70292 

58513 

46992 

36000 

25869 

16969 

9682 

4340 

1090 

0 

31200 

IO53OO 

58513 

46992 

54000 

19670 

25450 

9682 

4340 

1630 

0 

61833 

49977' 

38198 

26677 

15685 

5554 

-3346 

-00633 

-05975 

-09225 

-20315 

11.63 
u.63 . 

11.62 
11.61 
11.60 
11.58 
11.72 

11.63 
11.62 
11.62 
11.62 

1.000 

.808 

.617 

.431 

.254 

.090 

-.054 

-.172 

-.258 

-.311 

-.328 

356,777 


T 0 d) - -lil 


; x 341303 + 100 x 923* 3886000 


900 


200 


12 ' 


( 2 ) 


ffr x 356777 + 100 x 9682 4063000 


19, *30 


20,315 


® 2 # 11.62 (radians/ seo) 2 

KATIOKAL ADTTSORT 
COMOTTSE FOR ABtOM A PTT C S 


* 


« 


200 


200 


NACA TN No. 1522 



NACA. TN No. 1522 


55 


TABLE 11.- FIRST SYMMETRICAL BEIfDINO MODE OF A FREE-FREE BEAM 


CARRYING MASSES MOUNTED THROUGH 3PRINQS 
[l - 100 In. ; X - 10 In.; 1-5 in.*j E - 10,000,000 psi] 



1 

2 

3 

* 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 

15 

Sta- 

tion 

X 

m 

E 


E 

Tl (0) 

<») 

iff}/ 0 ) 

P eq 

V 


M 

Heq 

Slope 

Tj'd) 

Cosxo 

n 


125X 

n 

125*n 



_X_ 

12 

_x_ 

12 

12 

>£ 

12 

\3 

X3 

X* 

factors ' 

1** 

1** 

1**11 

1**311 

1***I 

10 

9 

8 

7 

6 

5 

* 

3 

s 

1 

0 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0. 38 
1.50 
1.00 
1.00 
1.50 

.76 

1.50 

1. ao 
1.00 
1.50 

.38 

1.000 

.855 

.712 

.571 

.*38 

.315 

.206 

.118 

.053 

.013 

0 

0.380 

1.282 

.712 

.571 

.658 

.239 

.309 

.118 

.053 

.019 

0 

0.769 

.62* 

.*81 

.3*0 

.207 

.08* 

-.025 

-.113 

-.178 

-.218 

-.231 

1.000 

.811 

.625 

.**2 

.269 

.109 

-.032 

-.1*7 

-.232 

-.283 

-.300 

5.62 

9.73 

7.50 

5.31 

3.2* 

1.33 

-.36 

-1.73 

-2.75 

-3.36 

-1.78 

*-02.65 

5.62 
15.35 
22.85 
28.16 
31. *0 
32.73 
32.37 
7.99 
5.2* 
1.88 

0 

5.6 

21.0 

*3.8 

72.0 

103.* 

136.1 

168.5 

176.5 
181.7 
IB3.6 

Tt 

259 

531 

867 

12*2 

1633 

9*6 

1029 

2115 

2177 

1098 

1197* 

11897 

11638 

11107 

102*0 

8998 

7365 

5390 

3275 

1098 

82982 

71008 

59111 

*7*73 

36366 

26126 

17128 

9763 

*373 

1098 

0 

11.52 

*.339 


‘ Y l (0) -Ti l(0) 

c 

vhere - — 


-T o ( 0) 

ff *( 0) + (-^'(oQsttttlon 3 

i -(sr 




0.231 


and p< 


, -5. _ 5000 , 
100 


50 ( radians/ sac) 2 and 


11 (assumed) 


16 

17 

IS 

19 

20 

Station 


TlCl) 

*i (0) A (1) 

Yjd) 

Common 

125X5. 

X* 



factors -5 '' 

1**®E1 

1**EI 



10 

31500 

63780 

11.28 

1.000 

9 

106600 

51810 

11.28 

.811 

8 

59111 

39910 

11.28 

.625 

7 

*7*73 

28270 

11.28 

.**3 

6 

5*500 

17170 

11.28 

.269 

5 

19630 

6926 

11.32 

.1083 

* 

25700 

-2072 

11.12 

-.032* 

3 

9763 

-9*37 

11.22 

-.1*78 

2 

*373 

-1*830 

11.25 

-.232 

1 

16*7 

-18100 

11.26 

-.28* 

0 

0 

360 .*97 

•-19200 

11.25 

-.301 


iS 5 i£ x 360*97 + 100 * 97 s 3 

T Cl) . -iiS 2i2« *382000 . 19200 

0 228.2 ® 2 8.2 v 

m 2 » 11.27 (radiaas/sec) 2 

For 1— Inch beam tip deflection, 

m^o 8 (Ti^Hstition 3 
Maximum spring force - i 

_ 100 x 11.28 x 0.1*78 . 2X6 it 

1 _ H- 88 

50 

HMI01WL AOTISOHT 
OOWITTKS FOR AER0HAWT1CS 





12,- jmat IBIXSIWWffilCAL HBHDIK More OP A PBn-BE* BUM OABHTma C0*CSaraA3ED MASSES 


CJi 

CT> 



* °- 017 ) + ( W >° *3 OX -O.6o) 

( 2 S? x 3<8V ) + [«»* (»)*] 


KABtOMl ADVISORY 

cot u r m pop aeromuttics 
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MBU 12.- FIRST AHTISYMNETHICAl BEHDIK3 RODE OF A FREE-FREE BEAK CARRYING CONCENTRATED MASSES - ConOlttJeO 
[l - 100 la. i X - 10 la.j 1-5 la.*; E - 10,000,000 pal] 


I 


* / 


18 

19 

20 

21 

22 

23 

24 

1 25 

26 

27 

28 

29 

30 

31 

38 

Station 

Slope 

gg| 

KXSl'f 1 ) 

£ 

r(l)x 

Y^l) 

EH 

P e9 

p aq 

3 

n 


Slope 

Yi'fS) 

Conaon v 

X3 

PS 

L 

12 5 XL 

■ 

wm 

■ 

_x_ 

12 

x_ 

12 

X_ 

12 

x£ 

12 

X3 

x3 

mm 


144EI 

144EI 

144 

144SI 

144EI 

144X1 

144SI 

I 

1 

3480 
3407 
3174 
2722 
. 2029 
, 1109 
0. 

^-1168 

-2027 

-2466 

10260 

6780 

3373 

199 

-2523 

-4552 

-5661 

-5661 

-4493 

-2466 

0 

10260 

6100 

2700 

139 

-3.515 

-2275 

-2265 

-1700 

-900 

-247 

0 

38^5 

9150 

2700 

139 

-2270 

-1730 

-3395 

-1700 

-900 

-370 

0 

5510 

-2880 

-2590 

-2300 

-2020 

-3.730 

-3.440 

-1150 

-860 

-570 

-290 

0 

13140 

9370 

5670 

2220 

-793 

-3110 

-4510 

-4800 

-3920 

-2170 

0 

1.000 

.713 

.431 

.169 

-.060 

-.236 

-.343 

-.365 

-.296 

-.165 

0 

5.42 

8.56 

5.19 

2.06 

^67 

-2.76 

-4.03 

-4.29 

-3.51 

-1.95 

0 

-43.80 

5.42 

13.98 

19.17 

21.23 

20.56 

17.80 

13.77 

—34.38 

-37.83 

-39.78 

0 

5.42 
19. 40 
38.57 
59.80 
80.36 
98.16 
111.93 
77.61 
39.78 
0 

74 

238 

465 

717 

962 

1174 

928 

475 

3630 

3556 

3318 

2853 

2136 

1174 

0 

-1231 

-8159 

-2634 

10643 

7013 

3457 

139 

-2714 

-4850 

-6024 

-6024 

-4793 

-2634 

0 


5510 125*6 + 100 x 30 x (-5661) 

s (1) _ lii 3- . S . aopOQ g -.28.8 

423000 X 23000 


33 

34 

35 

36 

37 

38 

39 

Station 


£ 

E< a )jc 

T^) 

VW) 

T^a) 

I Conaton \ Ti 

125XL 


x4 



factors 


144 


144SI 



10 

10643 

4045 

-3180 

13820 

58.1 

1.000 

9 

6310 

9460 

-2860 

9870 

51.9 

.714 

8 

2765 

2765 

-2540 

6000 

51.7 

.434 

7 

97 

97 

-2220 

2360 

51.5 

.171 

6 

-1629 

-2443 

-1910 

-804 

__ 

-.058 

5 

-2425 

-1843 

-1590 

-3260 

52.1 

-.236 

4 

-2410 

-3615 

-1270 

-4750 

52.0 

-.344 

3 . 

-1807 

-1807 

-950 

-5070 

51.9 

-.367 

2 

-959 

-959 

-640 

-4150 

51.7 

-.300 

1 

-263 

-394 

-320 

-2310 

51.4 

— .16T 

0 

0 

0 

0 

0 

— 

0 



5306 






'5306 i SS S t + 100 x 30 X (-6024) . 

s (2) _ ^ _ -13470000 _ 

423000 423000 


®1 2 t * 51.8 (radlaaa/aec) 2 
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WBU 13.- BKM WCn) OH 3JW30. FTB3T BEHDEH3 HOEQE 
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Figure 2,- Two successive numerical integrations of a curve by the summation method. 
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Figure 3 .- Formulas for equivalent concentrated loads. 
(Equations (a) and (b) from reference 5.) 
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Figure 5-" Satisfaction of boundary conditions for the first symmetrical mode of a free-free 
beam by the general method of approaoh. 
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Figure 6.- Satisfaction of boundary conditions for first antisymmetrical mode of a free-free 
beam by general method of approach. 
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Outline of Iteration Procedures 


(1) Deflection Y^ 0 ^ is assumed and loading mY^ 0 ^ is computed. 

(2) Shear is determined in two parts; known shear variation relative to left end, and an 

assumed correction for the left-end reaction. 

(3) By successive integrations, the moment, slope, and deflection diagrams are found 

from each of the two shear diagrams. 

( 4 ) Deflection y 0 must equal y' ; therefore, the deflection arising from the shear 

correction is adjusted. The two diagrams then added give the final deflection 
Y^ ' , which satisfies the deflection boundary conditions. 

(5) Process is repeated until Y^^ << Y^ + ^ . 


Figure 7 .- Boundary conditions and iteration procedure for a beam fixed at one 
end and simply supported at the other end. 
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Figure 8.- Graphical representation o 
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Figure 9.- Fundamental bending mode of a nonuniform 
cantilever beam. 



x/L 

Figure 10.- Fundamental torsional mode of a nonuniform 
cantilever beam. 
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Figure 12.- First symmetrical bending mode of a free-free 
beam carrying concentrated masses. 
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Figure 13.- Equations for correcting assumed deflection 
curve for a beam carrying spring- mounted masses. 
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Figure 1 I 4 ..- First anti symmetrical bending mode of a free-free 
beam carrying concentrated masses. 



Figure 15*~ Equation for correcting assumed 
deflection curve for a beam mounted on a 
spring . 


Relative deflection Relative deflection 


NACA TN No. 1522 


75 




Figure 16 . - First and second bending modes of a 
beam-mass-spring system. 



